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Production of x^-mesons at LHC 
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In the present paper we discuss the P-wave bottomonium production within 
both color octet and color singlet models in NLO at LHC energies. We calculate 
the cross section and the transverse momentum distributions for the Xi,2b states. 
We obtain, that the ratios of bottomonium and charmonium spin states production 
are fundamentally complementary at different p? scales. We give predictions for 
the ratio of the n = 2 and n = 1 radial excitations production cross sections. 

5^ PACS numbers.: 14.40.Pq, 14.40.Nd, 13.85.-t 

in : 1 Introduction 
> 

q^ ■ Studies of heavy quarkonia give a deeper understanding of the Quantum Chromodynamics. In 

the recent measurements made by ATLAS [1] and DO [2] Collaborations the new data for the bb 
systems was obtained. The radial excitation of the P-wave \b states was observed in radiative 
transitions to the 5*-wave T states, while Xfe(3P) was seen for the first time. The Xb system 
is a triplet of closely spaced states with total spin J = 0, 1 and 2, that were not observed 
separately at ATLAS yet. Usually, only 1 ++ and 2 ++ can be detected, since ++ state has very 
small radiative branching fraction. So, in the present paper, we shall focus on a production 

^ of 1 ++ and 2 ++ charmonium and bottomonium states. We shall consider the production with 

high transverse momentum of the final quarkonium. 

Main contribution to the production processes at high energies (~ TeV) is given by gluon- 
gluon subprocesses. The problem of the quarkonia production can be divided into two parts. 
The first part of the problem is to obtain the nonzero transverse momentum of the final 
quarkonia using the integrated partonic distributions (PDFs). The second part consists in 
the hadronization process, i.e. formation of the quarkonia with certain quantum numbers. For 
example, let us consider the production of the | 3 Si) (1 ) state. The colorless state with such 
quantum numbers cannot be produced from two gluons because of the charge parity conserva- 
tion. So, the production of the color singlet state with the additional emission of a gluon in 
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the final state was introduced [3, 4] to explain the quarkonia production. Additional emission 
of the gluon gives the p^-distribution and removes the prohibition of the 1 state production. 
It is well known, that such a model, in its naive understanding, is in the contradiction with the 
experimental data at high pt- In the case of P-wave meson production the situation is differ- 
ent: ++ and 2 ++ states can be produced from two gluons, but in the collinear approximation 
(with integrated over p? partonic distributions) the transverse momentum distribution of final 
quarkonium cannot be obtained. Moreover, | 3 Pi) state cannot be produced due to Landau- 
Yang theorem, which forbids the production of axial meson from two massless gluons. These 
problems have led to the approach, known as k<r- factorization [5, 6, 7, 8], in which unintegrated 
partonic distributions are used. However, we think that such an approach should play a signif- 
icant role at the low Pt region only, while in the high pt region the dominant contribution to 
the pt is given by the processes with emission of a single hard gluon in the final state. 

In this paper we try to resolve the above problems considering a set of diagrams, in which 
high pt is achieved by the emission of single hard gluon, mostly from the initial state. It turns 
out that all three P-waves states, including 1 ++ are allowed. As we shall see, the color-singlet 
contributions give pt dependence, which is similar to the experiment data for Xc mesons, for 
example from CDF [9], but the absolute normalization of the cross section is several times 
smaller than in the experiment. To resolve this contradiction we took into account the color- 
octet contributions to the cross sections. In the next section we show that if we consider the 
color-octet matrix elements as free parameters, the good agreement with the experimental data 
can be obtained for the case of Xc and the predictions for the Xb states can be provided. Next, 
we consider the cross section ratio of the 2 ++ and 1 ++ production both for Xc and Xb, which was 
recently measured for charmed mesons at the LHCb detector [10], and show, that our approach 
gives the value, lying at the lower limit of uncertainties of the LHCb data. There is no such 
experimental data on the bottomonium production yet and our result should be considered as 
the theoretical predictions for further experiments. Next, we consider the ratio of the n = 2 
and n = 1 radial excitations, which is not measured yet and present our theoretical predictions. 

2 Xc production at high pt 

The cross section of heavy quarkonium production in hadronic interaction can be expressed 
through the cross sections of hard subprocesses: 

da [pp ->• Q + X] = dx 1 dx 2 f g (x 1 ;Q) f g (x 2 ;Q)da[Q] (1) 

Here, Q is one of the | 3 Pj) states, x± y 2 are the momentum fractions of the partons, f g (2:1,2; Q) 
are the gluon distribution functions in initial hadrons, and da [Q] is the cross section of the 
meson production at the partonic level. 

As already mentioned in the Introduction, in order to obtain nonzero transverse momentum 
of final quarkonia using the usual partonic distributions, it is necessary to consider next to 
leading order subprocesses with the emission of an additional gluon, i.e. gg — > XbjQ- The 
corresponding Feynman diagrams are shown in Fig. 1. The partonic cross sections of these 





Figure 1: The Feynman diagrams of the gg — > Xb9 subprocesses. The first four diagrams are 
valid both for color-singlet and color-octet mechanism, while the last two diagrams corresponds 
to the octet production only. 



reactions were calculated by a number of authors [11, 12, 13, 14], but result presented in [13] 
disagrees with others works. The reason is that diagrams shown in Fig.l include 3-gluon vertex 
and the additional ghost contributions should be included. Our calculations performed in QCD 
axial gauge (which does not require additional ghosts contributions) reproduces [11, 12, 14] 
results. 

It should be noted, that the above expressions depend on two scale parameters: the renor- 
malization scale /i in strong coupling constant as(/i) and the maximum transverse momentum 
of collinear gluons Q in structure functions f(x,Q). At low energies the variation of these 
scale parameters leads to significant variation of the total cross sections [15], while for the 
energies about several TeVs the dependence is not crucial. In our work we will use the value 
Q = fi = M, where M is the mass of the heavy quarkonium. Also we shall neglect the mass 
difference between states from same xj{ n P) triplet. 

It is known, that color-singlet model of qq pair hadronization into observable meson does not 
fully describe experimental data. The reason is that color-singlet is only the first approximation 
in the Fock structure of quarkonium state [16]: 



\ 2S+1 Lj) = 0(l)\ 2S+1 LjW) 

+ 0(v)\ 2S+1 (L±l)j,Wg 

+ 0{v 2 )\ 2{ - s± V +1 Lj,Wg) 

+ 0(v 2 )\ 2S ^Lj^gg) 

+ ... 



(El) 
(Ml) 
(ElxEl) 



(2) 



where v is the relative velocity of quarks in heavy quarkonium. In the above expression El and 
Ml are electric (AL = 1, AS = 0) and magnetic (AS = 1, AL = 0) transitions respectively. 



In the case of \c mesons equation (2) gives: 

\Xcj) = ( 0^[ 3 Pf] ) \ 3 Pj [1] ) + ( 0**f sf] ) l 3 ^) + 

+ ( O^^P^} }\ 1 Pi [8] ) + J2( XcJ [ 3 Pf>] ) \ 3 P.J' [8] } + • • • (3) 

j> 

Here the fist two terms are order of 0(f), and next two terms are of order 0(v 2 ). 

In Fig. 2 the contributions from different states in (3) to the ^-distributions of J/i[) mesons 
produced via \c radiative decays: 

-r-\pP -> Xc + X -> J/rj) + X] = Bt[xci}-^-\pP -» Xa + X] + Br[x c2 \-^-\pp ->■ X c2 + A"]. 

are shown in comparison with the CDF experimental data [9]. The octet j 3 ^ ) state gives pt- 
distribution approximately described by ~ 1/py, which is far different from the experimental 
one, which, in turn, is approximately described by ~ 1/p^- So, in order to fit the experimental 
data, we should assume, that j 3 ^ ) gives a negligibly small contribution in the considered 
Pt interval and plays role for higher pt only. The other states from the Fock expansion (3), 
including color- singlet, have the p^-distribution similar to the experimental measurement. In 
order to fit the experimental data, we used CDF data for the Xc production and the available 
data for the cross section ratio of cr{x C 2) / v{x&) [10, 24, 25] (see next section and Fig. 4). The 
singlet contribution is determined by the singlet matrix element 0[ 3 Pj ], which is connected 
with the derivative of the radial part of the wave function at the origin: 

(O^[ 3 P [ J ] ]) = ^(2J+l)\R'(0)f. 

For the last one we took the value derived from the % c2 photonic width in LO: 

|i?'(0)| 2 = 0.075 GeV 5 . 
The octet matrix elements for \ 3 S[ ) satisfies the multiplicity relations 

(O x ^[ 3 Sf ] \) = (2J+l){O x ^[ 3 Sf ] }} 

For the last one we obtained from our fit the following value: 

(O Xc0 [ 3 Sf ] ]) < 2 x 10- 5 GeV 3 (4) 

Since ^Pi' 8 ]) and | 3 Pjt 8 l) states have similar p^-dependence, they could not be resolved 
separately and only a linear combination of the corresponding octet matrix elements can be 
determined. In the notations of [14] we obtain the following values for the non-perturbative 
parametres: 

(flXcofipWj) = (pxcoppjf)]) = o.oi GeV 5 (5) 
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Figure 2: Contributions to the Xc production from different states. Dotted line is color singlet 
contribution. Dashed line is a contribution from the octet | 3 PJ ) and ^Pj ; ) states. Dot- 
dashed line is color-octet \ 3 S[ ) state contribution. Solid line is a sum over all contributions. 
Experimental points are taken from CDF report [9]. 



The last result is independent on J'. Taking into account statistical weight of the Xcj-state the 
corresponding parametres for the Xcj production can be recovered: 

(^fpf 1 ]) = (2J+l)(P Xc0 [ 1 P 1 (8) ] 

(P*-fPJ? } ]) = (2J + l)(P Xc °[ 3 PJ?]) 

At this point we need to discusss the results of our fit. In most works (see, for example, 
review [17]) it is accepted to neglect the P-wave octet contributions due to the velocity power 
counting rules, and the following ranges of values of non-perturbative matrix element for the 
\ 3 S[ ) state are considered: 

(O Xc0 [ 3 Sf ] ]) ~ (0.1 - 0.5) x 10~ 2 GeV 3 

Nevertheless, we state that \ 3 S[ ) contribution gives inappropriate p^-dependence to describe 
the experimental data and therefore octet P-wave contributions should be included. 

On the other hand, the non-perturbative value of the octet \ 3 S[ ) contribution from work 
[8] is closer to our result: 

(O x <> [ 3 Sf ] })~2x 10~ 4 GeV 3 , 



but the derivative of the radial part of the wave function at the origin is about several times 
bigger then ours: 

|P'(0)| 2 ~0.37GeV 5 , 

so the almost entire contribution is given by the singlet term. We think, that this value is 
unreasonably high and does not agree with experimental data. Possible reason of the con- 
tradiction is that authors derive this value from the photonic width including QCD radiative 
correction. It is clear, that using such approach the radiative corrections to the initial state 
production should be considered on a par with radiative corrections in the final state, but this 
was not performed. 

3 Xb production at high £>t 

Let us proceed to the Xb mesons production. The picture of the bottomonium production is 
almost the same as in the charmonium case. The only difference is that probabilities in Fock 
expansion (3) have different values. We shall assume, that relative contributions from different 
Fock states are similar to charmonium mesons and do not depend on the radial excitation 
number n. So, for example, 



(0[ 3 P l f ] -> Xc ]) (Q[ 3 Pf ] -> XbinP)]) 
(OppW _> Xc] ) (oppW _> Xb{nP )]) 

From the dimension analysis, it is clear, that 

(Q[ 3 Sf ] ]) J_ 
{0[ 3 P [ j ] ]) ~ M2 

So, for \ 3 Si) we take 

(Q[ 3 Sf ] -> Xc]} _ Ml (OfS? -, Xb {nP)\) 
(0[*Pf -»- Xc}) M l (0[ 3 P? -»- Xb (nP)]) 



(6) 



(7) 



In Fig. 3 we show our results for the transverse momentum distributions for axial and tensor 
bottomonium. For the derivative of the wave function at the origin we take 

\R'J0)\ 2 = 1 GeV 5 , 



which is in agreement with potential models, presented in Tab. 1. For other non-perturbative 
parametres we use values obtained in the previous section and evaluated with (6) and (7). 
Presented curves allows one to easily reconstruct the distributions for any |i?' t (0)| value, by 
simple multiplying on the corresponding factor from Tab. 1. 

In Fig. 4 the p^-dependence of the Xq{J = ^)/Xq{J = 1) ratio is depicted. In case of 
charmonium mesons we compare our results with the LHCb data [10], CDF data [24] and CMS 
preliminary [25]. Both for (cc) and (66) mesons we present the cross section ratio in the rapidity 
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[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


Xb(lP) 


0.71 


1.02 


2.03 


0.94 


0.84 


0.51 


Xb(2P) 


0.99 


0.88 


2.2 


1.13 


0.98 


0.25 


Xb (3P) 


— 


— 


2.43 


1.17 


1.03 


— 



Table 1: \R' nP (0)| (in GeV 5 ) from different potential models 



> 

o 









lOOr 



lOr 



lr 



O.lr 



0.01- 



Xbi 2<y<4.5 

Xb2 2<y<4.5 

Xbi |y|<2 

Xb2 |y|<2 




Figure 3: Transverse momentum distributions of the Xbi and Xb2 mesons at s/s = 8 TeV. The 
derivative of the wave function at the origin is |i?^(0)| 2 = 1 GeV 5 , which is in agreement with 
potential models, presented in Tab. 1. Other non-perturbative parametres were obtained in the 
previous section in the case of charmonium and corresponds to the total curve in Fig. 2 and 
then evaluated for bottomonium using (6) and (7). 



range 2 < y < 4.5, but it is clear that this ratio depends weakly on the rapidity cut, so these 
predictions should be also valid for another rapidity regions. One can easily see, that within 
the experimental errors our results for charmonia agree well with experiment, so approach used 
in our paper is valid. It is clear, that since the cross section is proportional to the |i?'(0)| 2 
and if the mass gap between \J {nP) states with different radial numbers n is neglected, the 
ratio da[x2{nP)]/ da[xi{nP)\ does not depend on n and curves for charmonia and bottomonia 
presented on Fig. 4 are universal. 

One interesting property can be seen from Fig. 4: the bottomonium curve matches the 
charmonium curve if we perform the rescaling of the px variable: p? — > (M Xc /M Xb )pT ~ (1/3)pt- 
Moreover, this fact is exact model- independent theoretical result, which can be obtained from 
the dimension analysis. Indeed, the cross section of the quarkonia production depends on three 
dimensional parameters: hadronic energy, transverse momentum of final quarkonium and its 
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Figure 4: Transverse momentum distributions of the da [x2\ /da[xi] ratio. Solid and dashed 
lines stand for charmonium and bottomonium mesons. The dot-dashed line corresponds to 
the rescaled bottomonium ratio: abi/ab\{M Xc /M x ^pT)- As it is seen, it almost matches the 
charmonium curve. The experimental results for charmonium from LHC [10] are shown with 
dots, CDF [24] — with rectangles, and CMS [25] — with triangles. 



mass: 



da., 



da. 



(s,p T ,M) 



dpx dp? 

The cross section ratio J = 2 and J = 1 is dimensionless function of these parameters. Hence, 
we can write: 



da2(z 2 s,zpT,zM) j dai(z 2 s,zpx,zM) da 2 (s,pT, M) / dai(s,px,M) 



dp-i 



dp-i 



dpr I dpx 

Taking M = M Xc and z = M Xb /M Xc , we obtain 

da b2 (zp T ] s) I da bl (zp T ; s) da c2 {p T ] s/z 2 ) I da cl (p T ; s/z 2 ) 



(8) 



dpi 



dpi 



dpi 



dpi 



(9) 



Hence, if we know the charmonium ratio transverse momentum distribution at some cm. energy 
s, we can easily reconstruct the bottomonium ratio distribution at energy s (My/Mj/^) 2 . This 
relation is not violated by the scale dependence of a s and PDFs, since they are canceled in the 
ratio. 

It is intuitively clear, that at least for high energies the ratio weakly depends on s, since 
both axial and tensor meson acquire universal behaviour in s, so the s-dependence cancels in 
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this ratio. This argument does not work for low s, when another partonic channels becomes 
significant. The value of z is about z ~ 3, so the bottomonium distribution at y/s = 8 TeV 
corresponds to the charmonium distribution at y/s ~ 2.6 TeV, which is still very high, so the 
gluonic subprocess dominates. We can expect, that (9) can be approximately rewritten for high 
energies: 

da b2 {zp T ; s) I da bl (zp T ] s) da c2 (p T ] s) I da cl (p T \ s) 



dpx I dpx dpT I dpq 

3.1 Radial excitations 



(10) 



Let us now discuss the production of excited ^-mesons at LHC Recently Xb (lP)-, Xb (2P) — 
and %5(3P)-mesons were observed in T(lS)f and T(2S) , ~f modes by ATLAS Collaboration [1] 
and DO [2] Collaboration. In these experiments the radial excitation of the P-wave Xb states 
was observed in the radiative transitions to the S-wave T states. Unfortunately the exact values 
of the cross sections were not measured yet, but we can reconstruct them theoretically. 
The experimentally observable quantities are: 

2 

a th [nP,mS] = a th [pp -+ Xb (LP) + X -► T (mS) 7 + X] = J^ Br J [nP,mS] af [nP] (11) 

j=o 

where the following notations were used: 

af[nP] = a th [pp^ Xb j(nP)+X], 
BTj[nP,mS] = Br [ Xb j (nP) -)• T (IS) 7] . 

The branching fractions of radiative x&(lP)- and Xb(2P)-mesons decays are known experimen- 
tally [26]: 

B ri [IP, IS] = 35% ± 8%, Br 2 [IP, IS] = 22% ± 4%, 

Bn [2P, IS] = 8.5% ± 1.3%, Br 2 [2P, IS] = 7.1% ± 1%, 

Bn [2P, 2S] = 21% ± 4%, Br 2 [2P, 2S] = 16% ± 2.4%. 

In the case of scalar bottomonium the branching fractions are small, so we do not include it in 
the sum (11). 

In order to calculate the total cross sections, we should integrate the cross sections distri- 
butions obtained in the previous section and presented on Fig. 3 over appropriate pr region: 
A < p T < (s - M 2 )/2y/s. We use experimental values A = 12 GeV for ATLAS (\y\ < 2) [1], 
and A = 6 GeV for LHCb (2 < y < 4.5) [27]. Using these values we obtained the following 
prediction for the cross sections: 

I nn I Tin 

(LHCb) —-—a th [ Xbl (nP)] = 34.4-—, ^^ a th [ Xb2 (nP)] = 43- ' 



K(0)| 2 L "° 1V /J • GeV 5 \R' n (0)\ 2 GeV 5 

1 nn nn 

(ATLAS) —-—a^[ Xbl {nP)} = 5.2— . 7757777^^ ^2 (nP)] = 5.6- 



IK(0)| 2 L v n ' GeV 5 |P;(0)| 2 LA " V n ' GeV 5 



and for the ratio cr[x6(2P)]/o r [xb(lP)] in both cases: 



a 



th 



f^4l! = (0.29± 0.01 th ±0.1 6r ) 
IP, IS } K ' 



a th [IP, IS] 



lii 



21' 



R', 



w 



Here the first error is the error of our theoretical method, while the second is the uncertainty 
in the experimental values of branchings fractions. This ratio is determined by the ratio of 2P- 
and IP-states wave function derivatives at the origin (see Table 1). In Table 2 we show the 
predictions of different potential models for this ratio. 



Potential model 


Ratio Xb (2P)/ Xb (lP) 


[18] 


0.4 ± 0.01 th ± 0.14 6r 


[19] 


0.25 ± 0.01 th ± 0.1 6r 


[20] 


0.32 ± 0.02 th ± 0.1 6r 


[21] 


0.34 ± 0.01 th ± 0.12 fer 


[22] 


0.34 ± 0.01 th ± 0.12 fer 


[23] 


0.14 ± 0.05 6r 



Table 2: Theoretical predictions for the ratio of Xb{2P) and Xb(lP) production cross sections 
in different potential models 

In the case of x&(3P)-mesons production the corresponding branching fractions to T(LS') are 
not known yet, so we cannot perform similar analysis. However, we can estimate the radiative 
branching fraction of tensor meson using the following assumptions. We assume, that the main 
hadronic decay channel of tensor meson is 2-gluon decay: 



r[ X62 (3P) -V gg] 



128 o\R' 3 (0)\' 



-a 



5 ~° M 4 
and the radiative branching fraction of the x&2(3P) is equal to 

r[ X62 (3P) -► T(15) 7 ] 



Br 2 [3P, IS] 



T[Xb2(SP) -> T(l, 2, 35) 7 ] + r[x 62 (3P) -> gg] 



The radiative width can be obtained from the potential models [20, 22]. Using the cross section 
of T(LS') produced via x&i,2(3P) radiative decays 

a th [3P,lS] = Bn [3P, IS] af [3P] + Br 2 [3P, IS] erf [3P] , 

we can express the unknown branching fraction of the x&i(3P) — > T(1S')7 decay through the 

ratio of the 3P and IP production cross sections. The results are presented in Table 3, where 

we denoted 

_ a th [3P, IS] 

1 ~ a th [IP, IS] 

So, if the experimental value of 7 is found, the estimation on the radiative width x&l(3P) — y 
T(15')7 will be found too. Similar estimations can be also performed for the Xb -^ Y(2S)~f 
processes. 
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Potential model 


Br 2 [3P, IS] 


Bn [3P, 15] 


[20] 
[22] 


1.7% 
5.4% 


((41 ±9) 7 -1.3)% 
((42 ±9) 7 -4.1)% 



Table 3: Branching fractions of x&(3P) — > T(lS) , ~f decays 

4 Conclusions 

The present paper is devoted to the x^-meson production in hadronic experiments. Recently 
ATLAS [1] and DO [2] Collaborations reported about observation of Xb(l,2,3P) states, so 
theoretical description of these processes is desirable. 

In our paper we use both color singlet and color octet models for description of the P-wave 
heavy quarkonia production. For high energy reactions these mesons are produced mainly in 
gluon-gluon interactions and the cross sections of the corresponding processes can be written as 
a convolution of partonic cross sections and gluon distribution functions in initial hadrons. If 
one uses usual distribution functions integrated over the gluon transverse momentum and works 
at leading order (i.e. only gg — > Xb partonic reactions are considered), the information about 
the transverse momentum distribution of final quarkonia Pt is lost- It is clear, that in order to 
describe this distribution one has to use next to leading order results, when additional gluons 
in the final state are present. For high p? values processes with emission of additional gluons 
are suppressed by small strong coupling constant, so main contributions come from reactions 
with single hard gluon in the final state. For this reason we consider NLO partonic reactions 
gg —¥ XQJ9- It should be noted, that in this approach the production of axial Xi-meson is 
possible, while at leading order it is forbidden by Landau- Yang theorem. 

In our article we discuss transverse momentum distributions of Xc,b mesons. It is shown, that 
in the case of x c -meson production our predictions agree well with the available experimental 
data. In case of Xb{nP) mesons, we give predictions of their pT-distributions and absolute cross 
sections in ATLAS and LHCb experiments. One interesting property we found is that, accord- 
ing to our estimations, distributions of the ratios [da(x2)/dpT]/ [dcr{xi) / dpr\ for charmonium 
and bottomonium mesons coincide if the scaling p^ b — > (M Xc /M Xb )p^ b is performed. Using 
existing information about radial excitation of x^-mesons from different potential models we 
predict the ratio of x&(2P) and x&(lP) yields. Unfortunately, ATLAS and DO Collaborations do 
not provide the information about the normalization of their results, so currently comparison 
with experimental data is not possible. As for newly observed x&(3P) meson state, we show, 
that from the ratios of x&(3P) and x&(l,2P) production cross sections one can determine the 
branching fractions of x&i(3P) — > T(l, 2S)^ decays. 
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